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Abstract. We study solvable spin chains where either fields or couplings vary linearly 
in space and create a sandwich-like structure of the ground state. We find that the 
entanglement entropy between two halves of a chain varies logarithmically with the 
interface width. After quenching to a homogeneous critical system, the entropy grows 
logarithmically in time in the XX model, but quadratically in the transverse Ising 
chain. We explain this behaviour and indicate generalizations to other power laws. 



1. Introduction 

The entanglement between two parts of a quantum chain has been the topic of numerous 
recent studies [1]. For homogeneous systems, the entanglement entropy S in the ground 
state has been found to be a quantity of order one if the system is non-critical, while it 
varies as In L if it is critical. Here L is the length of the subsystem, which can be either 
an interval in an infinite chain or one-half of a finite chain. The constant in front of 
the logarithm is proportional to the number of contact points between the subsystems 
and to the central charge of the corresponding conformal field theory. In the vicinity 
of a phase transition, 5* varies as In^, where ^ is the correlation length, if L ^ ^ [2]. 
There have also been studies of non-homogeneous systems. For example, if a defect 
separates the two subsystems, the prefactor of InL varies with the defect strength in 
simple hopping models [3l H] , while for interacting electrons it scales either to zero or to 
the non- interacting value for large L [Sj E] . The central charge is also modified in chains 
with random couplings, but in this case by a constant factor of In 2 [3 El [9] . 

In the present study we consider non-homogeneous systems of a different nature. 
We assume that one parameter in the Hamiltonian varies linearly along the chain and 
consider two cases. In an XX model, we vary the strength of a magnetic field in the 
z-direction, while in a transverse Ising (TI) model we vary the couplings around the 
critical value. The first model constitutes a well-known problem, since it corresponds to 
free electrons hopping on a chain under the influence of a constant electric field. For a 
sufficiently large system, the central single-particle levels are then equidistant and form 
the famous Wannier-Stark ladder [ini [HI [El [iSl [H] . The situation has been realized 
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experimentally in optical lattices subject to a constant acceleration [T5l[TB] . The second 
model was investigated recently with regard to its critical properties [17]. Physically, 
they have in common that the gradient terms introduce an interface into the ground 
state, and a length scale A which measures its width. In the hopping model, this interface 
separates regions where the system is completely full and completely empty, respectively. 
In the TI model, it separates ordered and non-ordered regions. This interface should 
have a strong influence on the entanglement of the regions to the left and right of it. 
This is, indeed, what one finds. The entanglement entropy becomes constant for large 
L and the asymptotic value is determined by ln(A). Moreover, the deviation from the 
value in the homogeneous system has a scaling form in the variable (L/X). 

We also discuss the time evolution if the gradient is suddenly removed. For the 
hopping model, it turns out that S{\, t) depends only on the variable {t^ + A^) and 
thus can be obtained from the ground state entanglement with the field. This means 
in particular that it varies logarithmically in time. For the TI model, the situation is 
different and more interesting. Here one finds a quadratic increase of the entanglement 
with time. Such a behaviour has not been encountered before in such quenchs, but 
we show that it can be understood in a simple way using the quasiparticle picture of 
Calabrese and Cardy [T8] . 

In the following section 2 we review briefly the features of the Wannier-Stark 
problem. In section 3 we determine the correlation matrix from which S is calculated 
and present results for the XX chain in its ground state. In section 4 the TI model is 
considered, again in its ground state, while in section 5 results for the time evolution 
after the removal of the gradient are given for both models. In section 6 we sum up our 
findings and in an Appendix we comment on more general spatial inhomogeneities and 
the derivation of the length scales. 

2. Wannier-Stark problem 

The problem of lattice electrons in a homogeneous electric field has been the subject 
of many investigations. In the form of a simple one-dimensional tight-binding model it 
was studied, for example, in [191 [121 [131 [H] . The equivalent spin one-half XX chain was 
treated in [11]. The corresponding eigenvalue equation appears also in the treatment 
of certain reaction-diffusion models [20]. The Hamiltonian is, for a finite system of 2L 
sites with open ends, 

^ L-l L 

^ = ~2 ^ (4cn+i + cl+ic„) + /i {n-l/2)cicn, (1) 

n=-L+l n=-L+l 

Here the linear field is chosen such that it goes through zero between sites and 1 and 
thus is odd under a reflection of the chain. We will always assume h > 0. The eigenvalue 
equation for the single-particle states \k) then is 

- ^[(j)k{n - 1) + 0fc(n + 1)] + h{n - l/2)0fc(n) = uM^), (2) 
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and its general solution is given by a linear combination of the Bessel functions J„_k(1 / h) 
and Yn-n{l/h). The argument of these functions defines the characteristic length 
A = l/h, which will be of central importance in the following. For a finite system, 
K and ujk = h{K — 1/2) follow from the boundary conditions. The resulting spectrum is 
shown on the left of Fig. 1 for h = 0.05 and three values of L. One can see a linear region 
of equidistant levels with spacing h in the center, while the level separation becomes 
larger at the upper and lower end. Pictures of the corresponding eigenfunctions were 
first shown by Saitoh [12j. As L increases, the regions with nonlinear dispersion are 
moved towards ±oo and only the Wannier-Stark ladder with integer k = k remains. 
The eigenfunctions then are (j)k{n) = Jn-k{^/h) and concentrated near site k of the 
chain. In the ground state, the single-particle levels with A; < are occupied. The 
resulting density profiles are shown on the right of Fig. 1 for three values of the field. 
One can see that the transition from high to low density takes place in a region of width 
2A. 
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Figure 1. Wannier-Stark problem. Left: Single-particle eigenvalues uj^ for h = 0.05 
and three values of L. Right: Density profiles in the ground state for L = 50 and three 
values of the field. 



3. Entanglement in the XX chain 

In the following, we consider chains of 2L sites as in section 2 and study the entanglement 
between their left and right halves in the ground state. The corresponding entanglement 
entropy is = — Tr (p In p) where p denotes the reduced density matrix of one of the 
subsystems, e.g. the right half. Both p and S follow [211 [22] from the correlation matrix 

L 

Cmn = (clnCn) = ^ (l)k{m)(j)k{n)nk, (3) 

k=-L+l 

where the are the single-particle eigenfunctions appearing in ([2]) and the 

corresponding occupation numbers. In the ground state, these are one for A; < and 
zero otherwise. Restricting the matrix to the sites of the subsystem, 1 < m,n < L, and 
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calculating its eigenvalues (i, one obtains S as 

5 = - ^ InO - 5^(1 - 0) ln(l - Ci), (4) 

The calculation of the matrix and the diagonalization are done numerically. 

We first show that the length scale A introduced by the gradient appears directly 
in the single-particle eigenvectors of the correlation matrix. In Fig. [2] we have plotted, 
for fixed L, the eigenvector corresponding to the Q which is closest to 1/2 and thus 
gives the largest contribution to S. In the homogeneous case, it decays slowly from 
the boundary between the subsystems but extends through the whole interior. If the 
gradient is large enough, however, it becomes confined to a region near the boundary 
and the amplitude effectively vanishes at a distance A. This is similar to the situation 
in a homogeneous non-critical system, for example a dimerized hopping model [23]. In 
that case, the corresponding scale is the correlation length. 
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Figure 2. Eigenvector of the correlation matrix for L = 100 and five values of h. 
From top to bottom: h = 0, 0.01, 0.02, 0.05, 0.1. The center of the chain is at the left 
end. 

We now turn to the entanglement entropy. In Fig. [3] it is shown as a function of L 
for three values of h. For h = one has the well-known logarithmic increase [2] , but in a 
finite gradient S bends over after an initial rise and saturates rapidly. The change in the 
behaviour takes place if L ^ A. This saturation is easy to understand since the parts of 
the system outside the interface region are either full or empty and cannot contribute 
to the entanglement. 

The saturation value, calculated numerically for sizes L = 5A, is shown in Fig. [Has a 
function of A. It varies essentially logarithmically, with additional decaying oscillations. 
One can fit the data perfectly with the following form: 

5oo(A) = - ln(2A) + - + (5) 

where A ^ 1/4 and k = 0.726 is the non-universal constant appearing in the 
entanglement entropy of the homogeneous XX chain of length 2L, divided in the middle 

5hom = ^ln(4L/7r) + | (6) 
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Figure 3. Entanglement entropy for the XX chain with a linear field as a function of 
the half-length L 

The appearance of the constant in ([5]) is intriguing. One could make the two formulae 
identical by introducing an effective length Acs = n\/2 in (ISl), but this length would 
not have the simple interpretation of an interface width. In section 5.1 it will be seen 
that there is also a close relation of 6*00 (A) to the time- dependent entropy after a certain 
quench. 
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Figure 4. Asymptotic value 6*00 of the entanglement entropy as a function of the 
length A 

One should mention that in the case L —>■ 00 the correlation matrix can be written 
down analytically for any finite portion of the chain. As mentioned above, one then has 
to deal only with the Wannier-Stark ladder states, where the eigenf unctions are single 
Bessel functions of integer order. Thus 

00 

Cmn = Jk+mWJk+nW (7) 

k=0 





ln(2>.)/6 + 0.363 
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Using a sum rule, ([7]) can be rewritten as a simple product of Bessel functions. Dropping 
the arguments A it reads 

A 



Cri 



\Jm~lJn JmJi 



n-1 



In the limit A 



a 



l[m — n) 

oo, this reduces to the well-known result 
n)/2) 



sm vr m 



(8) 



(9) 



7[[m — n) 

of the homogeneous system. For general values of A it can easily be evaluated 
numerically. Since the Bessel functions become exponentially small for large order, 
taking L ^ A practically coincides with the result in the thermodynamical limit. 
Comparing with the numerical results using the exact eigenfunctions of ([1]) one also 
finds excellent agreement. 

We have also investigated the scaling limit, where the ratio of the length scales 
L/A = Lh is kept fixed while L ^ oo. We have calculated the entropy difference 
/S.S = S{L,h) — S{L,0) for several values of L. The results are shown in Fig. [5l 
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Figure 5. Scaling behaviour of the difference S{L, h) — S{L, 0) as a function of Lh for 
four system sizes. The full line in the right half is the function (|10p . 



In the case Lh ^ 1 we are in the Wannier-Stark limit discussed above, and for large 
L the curves tend towards a scaling function f{Lh) even for Lh values only slightly above 
1. This function is given by the difference of the asymptotic A ^ oo form of ([5]) and 
the formula ([6]) for the entropy of a half-chain. Taking the difference one finds 



6 2x 



(10) 



with X = Lh. The maximum is AS'max = /(I) = e f ~ 0.0753 in agreement with the 
value in the figure, and thus relatively small. It mirrors the slight bulge in the curve for 
S{L) seen in Fig. [3] before the asymptotic value is reached. 

For Lh < 1 the interface region fills the whole system, but the curves still show a 
nice scaling behaviour and AS" approaches zero quadratically as Lh 0. 
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Figure 6. Limiting value of the entanglement entropy in a TI chain as a function of 
X{g) — g^^^^. For comparison we have also plotted 1/12 Inx + ki with ki = 0.297. 

4. Transverse Ising chain 

We now consider the inhomogeneous quantum Ising chain with Hamiltonian 

L-l L 
n=-L+l n=~L+l 

where and are the components of a Pauli spin operator associated with site n, Jn 
is the nearest-neighbour exchange interaction and h the transverse field. The couplings 
are assumed to vary as 

J„ = J[1 + H- (12) 

The homogeneous chain with g = has a quantum critical point at J = h in the 
thermodynamic limit, L ^ oo. Thus setting h = J, the inhomogeneous system with 
g > is undercritical in the left half and overcritical, i.e. in the ordered phase, on the 
right. The length characterizing the transition region can be obtained from a scaling 
argument given in the Appendix [T7]. This leads to X{g) = ag"'^, for small g. Here 
uj = z//(l + z/) where u is the correlation length exponent. It enters because the 
perturbation drives the system away from criticality. With u = 1 for the TI model 
and choosing the constant a equal to one, \{g) = g~^^'^. The result for the XX chain 
can also be obtained in this way by using u = oo for the marginal perturbation one has 
there. 

The entanglement entropy between the two halves of the chain in its ground state 
is calculated again from the correlation functions by writing fllip in terms of fermions. 
Here one can work either in terms of Majorana operators [22j or of Fermi operators [9]. 
For S" as a function of L one finds the same overall behaviour as in the XX model, i.e. it 
rises logarithmically for small L and saturates then. The main difference is the absence 
of additional oscillations. This also holds for the asymptotic values which are plotted 
in Fig. [6] as a function of g~^^'^. The function Soo can be fitted by a simple logarithm 
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Figure 7. The entropy difference S(L,g)-S(L,0) as a function of a; = Lg^^^ for several 
values of L. The lines show the limiting behaviour, see text. 

5oo(A) = ^ln(2A) + A;i (13) 

The factor 1/12 instead of 1/6 corresponds to the central charge c = 1/2 of the TI 
model. The constant is given by ki = 0.297. In contrast to the XX case, there seems 
to be no relation to the constant 0.239 appearing in the entropy of the homogeneous 
model. 

With this information, we can again analyze the finite-size behavior of the entropy, 
which is expected to be 

S{L,g)-SmO) = f{L/X). (14) 

Indeed, for large L the entropy-difference approaches a universal function, as shown 
in FiglTl For small argument the scaling function behaves as S{x) ~ ~ g^, which 
follows from the fact that the entropy must be an even function of g. For large arguments 
the asymptotic behaviour is f{x) ~ — l/121na; + const which is in agreement with the 
result in Eq. (fT3ll . For intermediate values of x, the scaling function has a maximum at 
Xmax ~ 4.6 but, in contrast to the XX case in Fig. [5], is non-singular for all finite values 
of X. 

Instead of having the gradient in the couplings, one can also put it into the 
transverse field. In this case, due to duality, the ordered and disordered sides of the 
chain are reversed, but the size of the interface X{g) is expected from scaling theory to 
vary in the same way as before. Calculating the entanglement entropy, it is seen that 
S has again the form (fT3l) and also the scaling function in ( fT4l) approaches for large L 
a universal function with the same type of limiting behaviour. 

The three regimes of the system in a gradient, namely disordered (paramagnetic), 
interfacial and ordered (ferromagnetic) can also be probed by dividing the chain not in 
the middle but in two parts of length i and 2L — i and calculating the corresponding 
entanglement entropy. In Fig. [8] we present results for the case, where g = 1/L which 
means that the coupling J„ vanishes at the left boundary and equals 2 at the right one. 
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Figure 8. Entanglement entropy for TI chains divided asymetrically in two segments 
of length £ and 2L — £ for several system sizes and with g = 1/L. 



Both in the paramagnetic and in the ferromagnetic regime, the entropies approach 
a master curve as L is increased. In the interface regime, the entropy has a maximum 
which is of the same order as the value for i = L. The shape of the curves is reminiscent 
of those one finds in a homogeneous system of different finite sizes, if one varies the 
couphng constant [2]. In fact, if one uses the coupling Jhom(^) = — ^/L) in the 
homogeneous system and calculates 5* for a chain divided in the middle, the resulting 
curves are very close to those shown in the figure, even in the interface region, provided 
the size of the homogeneous system is close to the size of that region. This equivalence 
can be understood as follows. The length scale in the gradient system is A, thus the 
correlation between two points having a distance larger that A vanishes. Consequently, 
the contribution to the entropy at a given separation point x is coming from the sites in 
the range [x — A, x + A]. In this correlated domain, the couplings are varying only weakly 
and can be replaced by their average, which is just Jhom(^)- Therefore the entropy is 
that of the homogeneous system with size A and separation point in the middle. 

5. Time evolution after a quench 

The equilibrium results show that the main effect of the gradient is the introduction 
of a length scale, in terms of which the entropy still shows a logarithmic scaling. We 
now ask what happens if the gradient is suddenly switched off and a non-trivial time 
evolution of the state sets in. 

5.1. XX chain 

In the XX case, the fermionic operators evolve after the quench according to [21] 



where the sum is over the allowed momenta q = TTk/{2L + 1),A; = 1,2, ...2L for 
the homogeneous open chain, ipg{j) = L~^/^sin(g(j + L)) are the single-particle 




(15) 
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eigenfunctions and ujg = — cos q the corresponding eigenvalues. Therefore the correlation 
matrix C{t) at time t is obtained by multiplying C(0) from the left and right by the 
matrices U^(t) and U(t), respectively. The entropy is then calculated from C{t) as 
before. 

In Fig. [9] we show the resulting time evolution of S for various values of the initial 
gradient. For A = one is starting from a perfectly sharp domain wall. This situation 
has already been studied with regard to the evolution of the density in [25l |26] and 
with respect to the particle- number fluctuations in [27j. The entanglement entropy was 
obtained in a DMRG calculation in [28], but not investigated further. The curve looks 
very much like the one in Fig. HI on which we will comment presently. For larger A values, 
the initial entropy is higher, but the time evolution is also slower, and asymptotically 
all the curves seem to converge to the one with A = 0. Moreover, introducing the new 
variable r = \/W+X^ one finds an exact collapse of the data, as shown in the inset. 

1.4 
1.2 

1 

0.8 
^ 0.6 
0.4 
0.2 


20 40 60 80 100 

t 

Figure 9. Time evolution of the entanglement entropy in an XX chain with 
L = 150 after switching off the gradient. The initial values were, from top to bottom 
A = 50, 20, 10, 0. The inset shows S" as a function of the variable t = + 




This result can be derived analytically if one considers the limit L —> oo. In this 
case, one can work with a ring instead of an open chain and the quantities Ujm become 
Bessel functions. The equation for C(t) then reads explicitly 

Furthermore the matrix C(0) is given by ([7]) in this limit. Therefore one has sums 
of Bessel functions with two different arguments, t and A. Using their integral 
representations, the infinite sums over m and n can be carried out and one obtains 

CAt) = y^F;^,{X,t)F,^,{X,t) (17) 



fc>0 



where 



— 7T 
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In this integral one can now rewrite the variables as A = rcosy? and t = rsiny?. The 
addition theorem for the trigonometric functions then yields 

FniT,ip)= f ^e-'"^^"(«-'^)+*^" = e*^"J„(r). (19) 

with (f = arctan(t/A). 

Thus, up to phase factors the quantities F„ are Bessel functions with argument r, 
a fact which was not realized in [13]. At the level of the correlation matrix, the phase 
factors correspond to a simple unitary transformation C — U^CU and do not affect 
the entanglement entropy for which one finds the relation 

S{\,t) = 5(0, Vt2 + A2) (20) 

The entropy thus depends only on the variable r which is simply the distance from 
the domain-wall initial state in a space-time coordinate system where in the space 
direction we move to the equilibrium system with interface length A and then we 
further evolve this state in time. The lines of constant entropy are therefore circles 
in this quarter-plane. This explains in particular the result mentioned above that 
S'(A,0) = ^(O,^) for t —>■ X. Alternatively, for arbitrary time t one could think of the 
evolving state as being effectively the ground state of a gradient problem with interface 
length Aefr(t) = VX? + 1^- Hence, one has a front propagating with a time-dependent 
speed v{t) = dXcs/dt. 

Before closing this section we note that in finite systems one finds additional features 
at times r = 2L, 4L, . . ., which are larger than those in Fig. O Then the entropy shows 
a step-like increase which can be attributed to the propagating fronts which return to 
the center after being reflected at the open ends [2S]. 

5.2. TI chain 

In the TI chain, the calculation of the correlations using Majorana operators proceeds 
basically in the same way, but the time-dependent factors now contain the excitation 
energies Ug = 2sin(g/2) of the critical homogeneous system with h = J = 1. This gives 
a maximum velocity v = 1 for the excitations. In Fig. [10] the resulting entropy is shown 
for g = 1/1024 corresponding to A = 32, and four different lengths which were all much 
larger than A. One can see that 5* increases up to t ~ L and then drops again. The 
increase is different from the linear law one finds in homogeneous systems and can be 
fitted by a form 

S{t,g) = a{g)t'-b{g)t' + c{g) (21) 

where the cubic term represents a slight modification of the quadratic law. We have 
written a{g), b{g) and c{g), because for other g one finds the same time behaviour but 
with different coefficients. In particular a{g) is approximately linear in g, i.e. a{g) = ag 
where a ~ 0.23. 

The behaviour found above can be understood using the picture developed in [30] 
that at the quench pairs of quasiparticles are emitted which establish the entanglement 



Entanglement in spin chains with gradients 



12 



S(t) 




Figure 10. Time evolution of the entanglement entropy in the TI chain after switchin 
off a gradient g = 1/1024, for several values of L. The broken line is the fit using (|21[ 



between the parts of the system at later times. The number of these pairs must depend 
on the "distance" of the initial state from the ground state of the final system. This 
has to be connected with the difference between the initial coupling constants and the 
critical value J = 1. The precise expression can be found by considering a quench from 
a homogeneous initial system with J 7^ 1 and determining the coefficient a in the linear 
law S{t) — S{0) = at. This can be done numerically for our finite geometry, or by using 
the formula (2) in [31j for a segment in an infinite chain and dividing the result by 2. 
The latter approach gives a as an integral over momenta 



a(J) 



where Vr, 



and 



cos(g/2) is the velocity of the final quasiparticles, 
(J + 1) sin(g/2) 



V(J-l)' + 4Jsin2(g/2) 



H{y) = - 



y . A + y-, , 1 



y. 



2 ' 2 ' 2 '2 
Using partial integrations, this can be evaluated in closed form and gives 



a(J) 



\J-1\ 



TT 



— arctan( 



|J- 1| 
2y/j 



27rJ 



In 



J+1 



J -I 



(22) 



(23) 



(24) 



(25) 



For small | J — 1| this varies linearly, a ~ | J — l|/2, while for J ^ 00 (or for J ^ 0) 
it approaches the saturation value l/vr. The last term in (|25|1 is non-analytic at J = 1, 
but a good approximation in the range | J — 1| < 0.2 is 

a(J) = 0.4751 J- 1| -0.55(J- 1)2 (26) 

which we will use for simplicity. A similar result follows from the direct calculation in 
the finite geometry. 
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The quantity a/2 can now be used as the density of emitted pairs in the 
phenomenological formula of Calabrese and Cardy, which gives S{t), in a continuum 
approximation, as 



J -t 

where the integral counts all the pairs which end up in different halves of the system up 
to time t. In an inhomogeneous system, the obvious generalization of this formula is 



In the case = 1, this has exactly the form ( !2T1) and moreover the coefficient of the gt^ 
term equals 0.24 which is very close to the value 0.23 found by fitting. The agreement 
is also good for the cubic term where h{g)/g'^ equals 0.18 in (!29|) and 0.12 in the fit. 
Hence this formula describes the increase of S very well. In order to check it further, 
we have also studied the case 6 = 2, where the couplings increase quadratically from 
the center. Then S{t) looks very similar, but a closer analysis shows that it varies 
basically as gt^, which is again the prediction of (1291) . Also the coefficient has the 
correct value. Nevertheless one should point out that there is a region of small times 
where the numerical S{t) is rather fiat and not so well described by the formulae. This 
seems to hold up to t ~ A and would mean that inside the interface region the picture 
has to be modified. Another remark concerns the form of Eq. ( l28l) . To leading order, it 
is an integral over | J — 1 1 , which can also be interpreted as the (local) energy gap in the 
initial state. Such an expression was also used in a recent field theoretical treatment 
of inhomogeneous quenches, see Eq. (76) in [32] . In our approach we also obtain the 
exact prefactor on the lattice. However, the integral fl22|) from which it follows, does 
not permit to read off the result directly. 

Finally, let us comment on the decrease of S{t) beyond t = L. This is a feature 
of the finite geometry which one also finds in homogeneous quenches. There a linear 
increase of S is followed by an almost linear decrease and a zig-zag variation of S{t) 
results for larger times. It can also be understood in the quasiparticle picture. Due to 
the open ends, the quasiparticles moving outwards are refiected at the ends and follow 
their inward moving partners with a certain delay. As soon as they cross the middle, 
their contribution to the entanglement between left and right vanishes. This effect sets 
in when the quasiparticles from the ends arrive at the center, because their partners 
follow immediately. 

6. Summary 




(27) 




(28) 



If the initial couplings vary along the chain as 1 + gn , this gives 




(29) 



We have studied particular inhomogeneous systems where a power-law variation of some 
parameter introduces an interface with a certain width A. Beyond that region the ground 
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state approaches a product form. This suggests that instead of the real length of the 
system only A should enter the entanglement properties. In fact, we found that both 
for the XX model and the TI model the entropy in large systems is given by conformal 
expressions where A appears in the logarithms. This result is also plausible because one 
knows that the entanglement is connected with the interface between the two subsystems 
one considers. One finds the same In A behaviour in the g-symmetric XXZ Heisenberg 
chain, where the interface is produced by boundary fields [33] and the reduced density 
matrix is known explicitly [.34]. It is also somewhat similar to the situation for non- 
critical states, where the correlation length appears in the formulae. However, in our 
case, there is no translational invariance. For finite lengths L we have also shown that 
the entanglement entropy has a scaling form in the variable L/X. 

The time evolution after a removal of the gradient turned out to be particularly 
interesting. In the XX case, it is logarithmic as found usually in local quenches 
[211 [iHl [29] . Moreover, it displays a particular space-time symmetry relating static and 
dynamic entanglement. Formally, this results because both the single-particle states 
in the Wannier-Stark problem and the time evolution in the final homogeneous system 
are given by Bessel functions. Thus if one can treat, for example, the quench from a 
sharp interface, one has found the solution for all gradients and all times. However, one 
cannot apply the CFT approach of [35] since the walls there are of a different nature 
and lead to a linear behaviour of S{t). 

In the TI model, the time evolution after the quench turned out to be quadratic, a 
result not encountered before. We were able to explain this in the simple quasiparticle 
picture of [30] and could even give the numerical constants. We mentioned that more 
general power-law variations of the couplings lead to analogous results which can be 
understood in the same way. In a sense, this explanation works better than one might 
expect, because the assumption that quasiparticles are only emitted from nearby sites 
is not well founded near a critical point. Thus a more direct derivation following the 
lines of [32] would certainly be useful and interesting. 
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Appendix 

Although our main concern in this study were linearly varying parameters, similar results 
are obtained for other power laws. Consider a variation of the couplings in the transverse 
Ising model of the form 

Jn = J[l + 9\nf] . (30) 
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Figure 11. Entropy difference S{L,g) — S{L,0) as a function of Lg^^^ = L/\ for a 
quadratic variation of the couplings in the TI modeL 



This leads to a length scale A which can be estimated as follows. The typical deviation 
of the couplings from the critical value is A(A) ~ Jg^^, which leads to a length scale 
^ ~ A(A)^''. Since in the problem there is only one length scale, the width of the 
interface, we have ^ ~ A from which one obtains the self-consistency equation 

\r^[g\']-\ (31) 

with the solution: 

A = a . (32) 

Using V = 1 the exponent is = 1/(^^ + 1) which gives 1/2 for the linear variation 
considered in the main text. 

As an example, let us consider a quadratic variation, 6 = 2. Physically, this means 
that for J = h and g > the system is critical in the center and ordered more and 
more as one moves towards the ends. Thus one has a kind of sandwich structure with 
the thickness of the central part varying as A ~ g^^^^- Calculating the entropy, one 
finds again that it saturates for large L and the asymptotic value varies as in (fT3l) . 
The constant now has the value ^2 = 0.548 if one sets a = 1 in ( |32l) . The difference 
S{L,g) — S{L,0) is shown in Fig. [TTl The scaling function resembles that of the XX 
model in that it also has a cusp separating two different regimes. For small arguments 
it varies as x^, i.e. it is proportional to g. This behaviour is possible here since g < 
and g > are not equivalent. For large arguments it is given by — l/121na; + const. 

The case 6 = 1 corresponding to a linear "trapping potential" for the central part 
gives very similar results. The length is now A = g^^^ as in section 4, but the scaling 
function looks qualitatively as in Fig. [TTl However, for small arguments the behaviour 
is quadratic. 
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